Longitudinal single cell atlas identifies complex temporal relationship between type I interferon response and COVID-19 severity.

Due to the paucity of longitudinal molecular studies of COVID-19, particularly those covering the early stages of infection (Days 1-8 symptom onset), our understanding of host response over the disease course is limited. We perform longitudinal single cell RNA-seq on 286 blood samples from 108 age- and sex-matched COVID-19 patients, including 73 with early samples. We examine discrete cell subtypes and continuous cell states longitudinally, and we identify upregulation of type I IFN-stimulated genes (ISGs) as the predominant early signature of subsequent worsening of symptoms, which we validate in an independent cohort and corroborate by plasma markers. However, ISG expression is dynamic in progressors, spiking early and then rapidly receding to the level of severity-matched non-progressors. In contrast, cross-sectional analysis shows that ISG expression is deficient and IFN suppressors such as SOCS3 are upregulated in severe and critical COVID-19. We validate the latter in four independent cohorts, and SOCS3 inhibition reduces SARS-CoV-2 replication in vitro. In summary, we identify complexity in type I IFN response to COVID-19, as well as a potential avenue for host-directed therapy.

A single-cell atlas identifies pretreatment features of primary imatinib resistance in chronic myeloid leukemia.

Primary resistance to tyrosine kinase inhibitors (TKIs) is a significant barrier to optimal outcomes in chronic myeloid leukemia (CML), but factors contributing to response heterogeneity remain unclear. Using single-cell RNA (scRNA) sequencing, we identified 8 statistically significant features in pretreatment bone marrow, which correlated with either sensitivity (major molecular response or MMR) or extreme resistance to imatinib (eventual blast crisis [BC] transformation). Employing machine-learning, we identified leukemic stem cell (LSC) and natural killer (NK) cell gene expression profiles predicting imatinib response with >80% accuracy, including no false positives for predicting BC. A canonical erythroid-specifying (TAL1/KLF1/GATA1) regulon was a hallmark of LSCs from patients with MMR and was associated with erythroid progenitor [ERP] expansion in vivo (P < .05), and a 2- to 10-fold (6.3-fold in group A vs 1.09-fold in group C) erythroid over myeloid bias in vitro. Notably, ERPs demonstrated exquisite TKI sensitivity compared with myeloid progenitors (P < .001). These LSC features were lost with progressive resistance, and MYC- and IRF1-driven inflammatory regulons were evident in patients who progressed to transformation. Patients with MMR also exhibited a 56-fold expansion (P < .01) of a normally rare subset of hyperfunctional adaptive-like NK cells, which diminished with progressive resistance, whereas patients destined for BC accumulated inhibitory NKG2A+ NK cells favoring NK cell tolerance. Finally, we developed antibody panels to validate our scRNA-seq findings. These panels may be useful for prospective studies of primary resistance, and in assessing the contribution of predetermined vs acquired factors in TKI response heterogeneity.

Single-cell and bulk transcriptome sequencing identifies two epithelial tumor cell states and refines the consensus molecular classification of colorectal cancer.

The consensus molecular subtype (CMS) classification of colorectal cancer is based on bulk transcriptomics. The underlying epithelial cell diversity remains unclear. We analyzed 373,058 single-cell transcriptomes from 63 patients, focusing on 49,155 epithelial cells. We identified a pervasive genetic and transcriptomic dichotomy of malignant cells, based on distinct gene expression, DNA copy number and gene regulatory network. We recapitulated these subtypes in bulk transcriptomes from 3,614 patients. The two intrinsic subtypes, iCMS2 and iCMS3, refine CMS. iCMS3 comprises microsatellite unstable (MSI-H) cancers and one-third of microsatellite-stable (MSS) tumors. iCMS3 MSS cancers are transcriptomically more similar to MSI-H cancers than to other MSS cancers. CMS4 cancers had either iCMS2 or iCMS3 epithelium; the latter had the worst prognosis. We defined the intrinsic epithelial axis of colorectal cancer and propose a refined 'IMF' classification with five subtypes, combining intrinsic epithelial subtype (I), microsatellite instability status (M) and fibrosis (F).

Human Nasal Epithelial Cells Sustain Persistent SARS-CoV-2 Infection In Vitro, despite Eliciting a Prolonged Antiviral Response.

The dynamics of SARS-CoV-2 infection in COVID-19 patients are highly variable, with a subset of patients demonstrating prolonged virus shedding, which poses a significant challenge for disease management and transmission control. In this study, the long-term dynamics of SARS-CoV-2 infection were investigated using a human well-differentiated nasal epithelial cell (NEC) model of infection. NECs were observed to release SARS-CoV-2 virus onto the apical surface for up to 28 days postinfection (dpi), further corroborated by viral antigen staining. Single-cell transcriptome sequencing (sc-seq) was utilized to explore the host response from infected NECs after short-term (3-dpi) and long-term (28-dpi) infection. We identified a unique population of cells harboring high viral loads present at both 3 and 28 dpi, characterized by expression of cell stress-related genes DDIT3 and ATF3 and enriched for genes involved in tumor necrosis factor alpha (TNF-α) signaling and apoptosis. Remarkably, this sc-seq analysis revealed an antiviral gene signature within all NEC cell types even at 28 dpi. We demonstrate increased replication of basal cells, absence of widespread cell death within the epithelial monolayer, and the ability of SARS-CoV-2 to replicate despite a continuous interferon response as factors likely contributing to SARS-CoV-2 persistence. This study provides a model system for development of therapeutics aimed at improving viral clearance in immunocompromised patients and implies a crucial role for immune cells in mediating viral clearance from infected epithelia. IMPORTANCE Increasing medical attention has been drawn to the persistence of symptoms (long-COVID syndrome) or live virus shedding from subsets of COVID-19 patients weeks to months after the initial onset of symptoms. In vitro approaches to model viral or symptom persistence are needed to fully dissect the complex and likely varied mechanisms underlying these clinical observations. We show that in vitro differentiated human NECs are persistently infected with SARS-CoV-2 for up to 28 dpi. This viral replication occurred despite the presence of an antiviral gene signature across all NEC cell types even at 28 dpi. This indicates that epithelial cell intrinsic antiviral responses are insufficient for the clearance of SARS-CoV-2, implying an essential role for tissue-resident and infiltrating immune cells for eventual viral clearance from infected airway tissue in COVID-19 patients.

RCA2: a scalable supervised clustering algorithm that reduces batch effects in scRNA-seq data.

The transcriptomic diversity of cell types in the human body can be analysed in unprecedented detail using single cell (SC) technologies. Unsupervised clustering of SC transcriptomes, which is the default technique for defining cell types, is prone to group cells by technical, rather than biological, variation. Compared to de-novo (unsupervised) clustering, we demonstrate using multiple benchmarks that supervised clustering, which uses reference transcriptomes as a guide, is robust to batch effects and data quality artifacts. Here, we present RCA2, the first algorithm to combine reference projection (batch effect robustness) with graph-based clustering (scalability). In addition, RCA2 provides a user-friendly framework incorporating multiple commonly used downstream analysis modules. RCA2 also provides new reference panels for human and mouse and supports generation of custom panels. Furthermore, RCA2 facilitates cell type-specific QC, which is essential for accurate clustering of data from heterogeneous tissues. We demonstrate the advantages of RCA2 on SC data from human bone marrow, healthy PBMCs and PBMCs from COVID-19 patients. Scalable supervised clustering methods such as RCA2 will facilitate unified analysis of cohort-scale SC datasets.

Spearmint R2R3-MYB transcription factor MsMYB negatively regulates monoterpene production and suppresses the expression of geranyl diphosphate synthase large subunit (MsGPPS.LSU).

Many aromatic plants, such as spearmint, produce valuable essential oils in specialized structures called peltate glandular trichomes (PGTs). Understanding the regulatory mechanisms behind the production of these important secondary metabolites will help design new approaches to engineer them. Here, we identified a PGT-specific R2R3-MYB gene, MsMYB, from comparative RNA-Seq data of spearmint and functionally characterized it. Analysis of MsMYB-RNAi transgenic lines showed increased levels of monoterpenes, and MsMYB-overexpressing lines exhibited decreased levels of monoterpenes. These results suggest that MsMYB is a novel negative regulator of monoterpene biosynthesis. Ectopic expression of MsMYB, in sweet basil and tobacco, perturbed sesquiterpene- and diterpene-derived metabolite production. In addition, we found that MsMYB binds to cis-elements of MsGPPS.LSU and suppresses its expression. Phylogenetic analysis placed MsMYB in subgroup 7 of R2R3-MYBs whose members govern phenylpropanoid pathway and are regulated by miR858. Analysis of transgenic lines showed that MsMYB is more specific to terpene biosynthesis as it did not affect metabolites derived from phenylpropanoid pathway. Further, our results indicate that MsMYB is probably not regulated by miR858, like other members of subgroup 7.

Metabolic engineering of terpene biosynthesis in plants using a trichome-specific transcription factor MsYABBY5 from spearmint (Mentha spicata).

In many aromatic plants including spearmint (Mentha spicata), the sites of secondary metabolite production are tiny specialized structures called peltate glandular trichomes (PGT). Having high commercial values, these secondary metabolites are exploited largely as flavours, fragrances and pharmaceuticals. But, knowledge about transcription factors (TFs) that regulate secondary metabolism in PGT remains elusive. Understanding the role of TFs in secondary metabolism pathway will aid in metabolic engineering for increased yield of secondary metabolites and also the development of new production techniques for valuable metabolites. Here, we isolated and functionally characterized a novel MsYABBY5 gene that is preferentially expressed in PGT of spearmint. We generated transgenic plants in which MsYABBY5 was either overexpressed or silenced using RNA interference (RNAi). Analysis of the transgenic lines showed that the reduced expression of MsYABBY5 led to increased levels of terpenes and that overexpression decreased terpene levels. Additionally, ectopic expression of MsYABBY5 in Ocimum basilicum and Nicotiana sylvestris decreased secondary metabolite production in them, suggesting that the encoded transcription factor is probably a repressor of secondary metabolism.

A terpenoid phytoalexin plays a role in basal defense of Nicotiana benthamiana against Potato virus X.

Terpenoid phytoalexins function as defense compound against a broad spectrum of pathogens and pests in the plant kingdom. However, the role of phytoalexin in antiviral defense is still elusive. In this study, we identified the biosynthesis pathway of a sesquiterpenoid phytoalexin, capsidiol 3-acetate as an antiviral response against RNA virus Potato Virus X (PVX) in Nicotiana benthamiana. NbTPS1 and NbEAH genes were found strongly induced by PVX-infection. Enzymatic activity and genetic evidence indicated that both genes were involved in the PVX-induced biosynthesis of capsidiol 3-acetate. NbTPS1- or NbEAH-silenced plant was more susceptible to PVX. The accumulation of capsidiol 3-acetate in PVX-infected plant was partially regulated by jasmonic acid signaling receptor COI1. These findings provide an insight into a novel mechanism of how plant uses the basal arsenal machinery to mount a fight against virus attack even in susceptible species.

Over-expression of OSRIP18 increases drought and salt tolerance in transgenic rice plants.

Both drought and high salinity stresses are major abiotic factors that limit the yield of agricultural crops. Transgenic techniques have been regarded as effective ways to improve crops in their tolerance to these abiotic stresses. Functional characterization of genes is the prerequisite to identify candidates for such improvement. Here, we have investigated the biological functions of an Oryza sativa Ribosome-inactivating protein gene 18 (OSRIP18) by ectopically expressing this gene under the control of CaMV 35S promoter in the rice genome. We have generated 11 independent transgenic rice plants and all of them showed significantly increased tolerance to drought and high salinity stresses. Global gene expression changes by Microarray analysis showed that more than 100 probe sets were detected with up-regulated expression abundance while signals from only three probe sets were down-regulated after over-expression of OSRIP18. Most of them were not regulated by drought or high salinity stresses. Our data suggested that the increased tolerance to these abiotic stresses in transgenic plants might be due to up-regulation of some stress-dependent/independent genes and OSRIP18 may be potentially useful in further improving plant tolerance to various abiotic stresses by over-expression.

Genome-wide survey of the RIP domain family in Oryza sativa and their expression profiles under various abiotic and biotic stresses.

Ribosome-inactivating proteins (RIPs) are N-glycosidases that inhibit protein synthesis by depurinating rRNA. Despite their identification more than 25 years ago, little is known about their biological functions. Here, we report a genome-wide identification of the RIP family in rice based on the complete genome sequence analysis. Our data show that rice genome encodes at least 31 members of this family and they all belong to type 1 RIP genes. This family might have evolved in parallel to species evolution and genome-wide duplications represent the major mechanism for this family expansion. Subsequently, we analyzed their expression under biotic (bacteria and fungus infection), abiotic (cold, drought and salinity) and the phytohormone ABA treatment. These data showed that some members of this family were expressed in various tissues with differentiated expression abundances whereas several members showed no expression under normal growth conditions or various environmental stresses. On the other hand, the expression of many RIP members was regulated by various abiotic and biotic stresses. All these data suggested that specific members of the RIP family in rice might play important roles in biotic and abiotic stress-related biological processes and function as a regulator of various environmental cues and hormone signaling. They may be potentially useful in improving plant tolerance to various abiotic and biotic stresses by over-expressing or suppressing these genes.

A comprehensive transcriptional profiling of the WRKY gene family in rice under various abiotic and phytohormone treatments.

WRKY transcription factors play important roles in the regulation of various biological processes. We have analyzed the publicly available rice genome sequence databases and predicted 103 genes encoding WRKY transcription factors. Among them, the majority of rice WRKY genes (77.7%) were located in duplicated regions; 45.6% of WRKY genes were fragmentally duplicated and 35% of them were tandemly duplicated. These results suggested that genome duplications might be regarded as a major mechanism for expansion of this family in the rice genome. Subsequently, we analyzed their expression profiles under normal and abiotic stress, as well as various hormone treatments. Under normal growth conditions, 65 WRKY genes were expressed differentially either in their transcript abundance or in their expression patterns. Under abiotic (cold, drought and salinity) stresses and various phytohormone treatments, 54 WRKY genes exhibited significant differences in their transcript abundance; among them three genes were expressed only in stressed conditions. Among the stress-inducible genes, 13 genes were regulated only by abiotic stresses, another set of 13 genes were responsive to only phytohormone treatments and the remaining 28 genes were regulated by both factors, suggesting an interaction between abiotic stress and hormone signaling. On the other hand, we have also surveyed the expression divergence of duplicated genes under normal or stressed conditions, and the results showed that high expression divergence has occurred not only among fragmentally but also among tandemly duplicated genes. These results suggested that the high expression divergence could be one of the mechanisms for the retention of these duplicated WRKY genes.

Ds insertion mutagenesis as an efficient tool to produce diverse variations for rice breeding.

The availability of diversified germplasm resources is the most important for developing improved rice varieties with higher seed yield or tolerance to various biotic or abiotic stresses. Here we report an efficient tool to create increased variations in rice by maize Ac/Ds transposon (a gene trap system) insertion mutagenesis. We have generated around 20,000 Ds insertion rice lines of which majority are homozygous for Ds element. We subjected these lines to phenotypic and abiotic stress screens and evaluated these lines with respect to their seed yields and other agronomic traits as well as their tolerance to drought, salinity and cold. Based on this evaluation, we observed that random Ds insertions into rice genome have led to diverse variations including a range of morphological and conditional phenotypes. Such differences in phenotype among these lines were accompanied by differential gene expression revealed by GUS histochemical staining of gene trapped lines. Among the various phenotypes identified, some Ds lines showed significantly higher grain yield compared to wild-type plants under normal growth conditions indicating that rice could be improved in grain yield by disrupting certain endogenous genes. In addition, several 1,000s of Ds lines were subjected to abiotic stresses to identify conditional mutants. Subsequent to these screens, over 800 lines responsive to drought, salinity or cold stress were obtained, suggesting that rice has the genetic potential to survive under abiotic stresses when appropriate endogenous genes were suppressed. The mutant lines that have higher seed yielding potential or display higher tolerance to abiotic stresses may be used for rice breeding by conventional backcrossing combining with molecular marker-assisted selection. In addition, by exploiting the behavior of Ds to leave footprints upon remobilization, we have shown an alternative strategy to develop new rice varieties without foreign DNA sequences in their genome.

Genome-wide analysis of cyclin family in rice (Oryza Sativa L.).

The cyclins together with highly conserved cyclin-dependent kinases regulate cell cycle progression in plants. Although extensive and systematic study on cell cycle mechanisms and cyclin functions in yeasts and animals has been carried out, only a small number of plant cyclins have been characterized and classified functionally and phylogenetically. We identified several types of cyclin genes in the rice genome and characterized them by phylogenetic, tandem and segmental duplications analyses. Our results indicated that there were at least 49 predicted rice cyclin genes in the rice genome, and they were distributed on 12 chromosomes. Of these cyclins, one possessed only cyclin_C domain and no cyclin_N domain, and the remaining 48 cyclins with cyclin_N domains were classified as nine types based on evolutionary relationships. Eight of these nine types were common between rice and Arabidopsis, whereas only one, known as F-type cyclins, was unique to rice. No homologues of the F-type cyclins in plants could be retrieved from the public databases, and reverse transcription-PCR analysis supported an existence of the F-type cyclin genes. Sequence alignment suggested that the cyclin genes in the rice genome experienced a mass of gene tandem and segmental duplications occurred on seven chromosomes related to the origins of new cyclin genes. Our study provided an opportunity to facilitate assessment and classification of new members, serving as a guide for further functional elucidation of rice cyclins.